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ABSTRACT 

We develop a Lagrangian perturbation framework for rotating non-relativistic superfluid neu- 
tron stars. This leads to the first generalization of classic work on the stability properties of 
rotating stars to models which account for the presence of potentially weakly coupled su- 
perfluid components. Our analysis is based on the standard two-fluid model expected to be 
relevant for the conditions that prevail in the outer core of mature neutron stars. We discuss 
the implications of our results for dynamical and secular instabilities of a simple neutron star 
model in which the two fluids are allowed to assume different (uniform) rotation rates. 
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1 INTRODUCTION 

' Investigations into the stability properties of rotating self- 
gravitating bodies are of obvious relevance to astrophysics. By im- 
proving our understanding of the relevant issues we can hope to 
shed light on the nature of the various dynamical and secular in- 
stabilities that may govern the spin-evolution of rotating stars. The 
relevance of such knowledge for neutron star astrophysics may be 
highly significant. It is for example not clear (i) whether instabili- 
. ties may affect nascent neutron stars leading to the slow birth spin 
' rates inferred for radio pulsars, (ii) to what extent gravitational- 
wave driven instabilities may provide the explanation for the ob- 
served clustering of spin periods in Low-Mass X-ray Binaries, and 
(iii) if compact stars are likely to evolve through phases where un- 
stable oscillations lead to the emission of a detectable gra vitational- 
wave signal (see the review article bv lAnderssonI i2003h for a dis- 
cussion of these issues and further references). 

The aim of this (and a subsequent) paper is to develop a La- 
grangian perturbation framework for rotating non-relativistic super- 
fluid neutron stars. This leads to the first generalization of Friedman 
and Schutz's classic work ( Friedman & Schutz 1978a b) on the sta- 
bility properties of rotating stars to models which account for the 
presence of potentially weakly coupled superfluid components. We 
believe that our extension of the, by now 25 years old, single fluid 
results corresponds to a significant step forwards for this area of 
research. 

Ou r analy si s is based on the standa rd two-fluid model 
jCarteil Il989|- ICo mer & Lanelois' "1994"; "Carter & Lanelois 
_1998t iLanglois e t al. 1198; Lindblom & Mendefl 2000 ; 
Sedrakian & Wasserman _|^Kdt Andersson & Comer 120011; 
Comeh.20 02; Lee & Yoshidijl |2003l lYoshida & Lee..2003 JbtlPrixl 
2004 ICart er & Chamel 20031® expected to be relevant for the 
conditions that prevail in the outer core of mature neutron stars. In 



order to keep the analysis tractable we are forced to make many 
simplifying assumptions, but our model nevertheless includes the 
extra degrees of freedom associated with any two-fluid config- 
uration. Most importantly, we allow our stationary background 
configurations to be such that the two fluids rotate at different 
rates. This is a key addition to the standard picture. It is a necessary 
inclusion since all seemingly viable models of the large Vela pulsar 
glitches are based on the notion of angular momentum transfer 
between different fluid components rotating at different rates. It 
is also an exciting addition to the standard model, because the 
additional degree of freedom may lead to new phenomena. An 
interesting example of this is the recently discovered superfluid 
two-stream instability fAndersson et aL .2003^ .Prix et aL .20041 
lAndersson et all2002l() . 

The layout of this paper is as follows: We begin by introduc- 
ing the two-fluid model for superfluid stars and the Lagrangian 
perturbation theory that we will employ (Section |2j. We then re- 
visit (in Section |3j what should be well-known results for a sin- 
gle (barotropic) fluid model. We also revise those results in such a 
way that the main concepts will remain relevant for the two-fluid 
problem. Yet they are easier to introduce in the single fluid case. 
Having developed the necessary tools we then move on to the two- 
fluid case (Section |4}- As we will see, the extension of the single 
fluid formulae is relatively straightforward provided that we only 
allow for the direct chemical coupling betwee n the two fluids [for 
the m oment neglecting the entrainment effect jAnd reev & Bashkinl 
Il975t Borumand et al ]|l996tEomer & Jovntll2003HComedl2004b . 
the role of which will be the main focus of a subsequent paper]. 
Section |5| then provides a discussion of some implications of our 
results, which is followed in Section|6|by specific implications for 
the r-mode instability and gravitational wave emission. Finally, in 
SectionQwe briefly look ahead to the work that will be presented 
in a subsequent paper. 
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2 TWO KEY INGREDIENTS 

We begin our discussion by introducing the two main ingredients 
of our analysis: The two-fluid model for ne utron stars, and the 
Lagrangian perturbation theory developed bv iFriedman & Schutd 
a978a.b.) . 



2.1 The two-fluid model for superfluid stars 

A mature neutron star is likely to contain several superfluid com- 
ponents. One would certainly expect superfluid neutrons to coexist 
with the crystal lattice of nuclei that makes up the outer kilometer or 
so of the star. In the outer parts of the fluid core superfluid neutrons 
are expected to coexist with superconducting protons, while the 
deep core may contain exotic states of matter like superfluid hyper- 
ons and perhaps even colour-superconducting deconfined quarks. 
The modelling of a star comprising all these components is far be- 
yond our current means. In fact, the parameters governing the pos- 
sible states are uncertain to say the least. Given this, most investiga- 
tions into the dynamics of superfluid neutron stars have considered 
a much simplified scenario that accounts for the presence of two 
distinct fluid components. 

The two-fluid picture of superfluid neutron stars is based on 
the notion that the outer core is dominated by superfluid neutrons, 
superconducting protons and electrons. Since the charged compo- 
nents couple electromagnetically on a short timescale they are as- 
sumed to move together. This assumption should be valid as long 
as we consider dynamics that takes place on a timescale long com- 
pared to that associated with the electromagnetic interaction. The 
outcome is a model describing the motion of the neutrons and the 
"protons" (a conglomerate of the charged components). 

The equations that describe the two-fluid model are derived 
from an energy functional E{nn,np,w^) where Wx are the two 
number densities. Given the two fluid velocities vf the quantity 

= — vf (Wyx = Wy — ui) represents the relative velocity be- 
tween the two fluids, for ease of notation we define lu^ = Wy^w'^^. 
Throughout this paper we will use the constituent indices x = n, p 
and y 7^ X. This allows us to condense the various equations sig- 
nificantly, but in order to avoid confusion we should stress that re- 
peated constituent indices never imply summation, while repeated 
vector component indices always do. It is also worth pointing out 
that all calculations will be carried out in a coordinate basis. Hence, 
we will distinguish between co- and contravariant vectors etcetera. 

Once we are supplied with the energy functional (the equation 
of state) a straightforward variation yields 



dE — tJ-xdn^^ + adw^ 

X— n.p 



(1) 



This leads to the definition of the two chemical potentials 
Mx^(l^) . (2) 



For later convenience, we introduce the notation /ix = ^^/m^, 
where we take the two nucleon masses to be equal such that rriB = 
rrin ~ nip. We also see that the entrainment coefficient a follows 
from 



( — ) 



(3) 



The dynamical equations that govern the two coupled fluids 
can b e derived either from a Newto nian variational p rinciple I Prix 
|2004) or from the Newtonian limit i Andersson & Comeil200L') of 



the f ully relat ivistic equations der ived by Carter and collab orators 
lCa rteJll989tJCo mer & Laneloi3 fT99i ICarter & Laneloi j fT993: 
iLanglois et alll998li) . The end result is two continuity equations 

9tnx + Vi(nx«j) =0 (4) 

and two coupled Euler equations 

(dt + viV,) (vt + e^wD + V. ($ + ^x) 

+£xiff V,ui = . (5) 

Here we have defined £x ~ 2a /n^. We also have the standard 
Poisson equation for the gravitational potential $; 



V $ = 47rmBG (Tin + Up) 



(6) 



We assume the background to be stationary and axisymmetric, 
with the two fluids rotating around the z-axis with rates Qn and fip 
respectively. Hence we have 



Ux ~ f^xV^' and lUyx = (^y ^ ^^x)v'' 
with ip^ given by 



(7) 



(8) 



In spherical coordinates, i.e. £ {r, 9, ip} , this vector has the 
components ip^ — (0, 0, 1), and its norm is tp'ipi = sin^ 9. 

In order to simplify the analysis, we restrict our attention to 
models with uniform rotation, i.e. take S7x to be constants. Further- 
more, in order to elucidate the details we will only consider the case 
of vanishing entrainment. That is, we let a — 0. The motivation for 
this is simple: If we could not make progress even in this restricted 
case, it would be pointless to consider the much more complicated 
problem which includes the entrainment. However, as we will see, 
the case of vanishing entrainment works out very neatly. Hence 
there is every reason for optimism, and we will return to the case 
Q 7^ in a subsequent paper. 

2.2 Lagrangian perturbation theory 

Our aim is to derive conserved quantities for the perturbations of 
the system of equations described above. By doing this we hope to 
be able to provide criteria that can be used to decide when the oscil- 
lations of a rotating superfluid neutron star are unstable. The need 
for such criteria is clear given that i) the astrophysical relevance 
of the gravitational-wave driven r-mode instability may to a large 
extent depend on whether the star contains superfluid components 
(Andersson 2003), and ii) the dynamical two -stream instability may 
set in above a critical relative rota tion rate I Andersson et al. 2003( 
IPrix et all200llAndersson et all2002l) . A detailed analysis of both 
these problems clearly requires an improved understanding of the 
stability properties of superfluid stars. 

Following Friedman & Schutz 1 1978a b), we analyze the 
problem within the Lagrangian perturbation formalism. The La- 
grangian variation AQ of a quantity Q is related to the Eulerian 
variation SQ hy 



AQ = SQ+£^Q 

where the Lie derivative £^ has the meaning 
for scalars, 

for contravariant vectors, and 



(9) 



(10) 



(11) 
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£e,Vi^(,^VjVi + VjSJ^i^ (12) 

for covariant vectors. 

The Lagrangian change in the fluid velocity follows from 

A^;' = dtC (13) 

where 5' is the Lagrangian displacement. Given this, and 

Agzj = V,Cj + VjC> (14) 
we have 

/^vi = dtii + v^Viij + -"^^i^i ■ (15) 
It is also useful to note that the Eulerian variations are given 



by 
and 



(16) 



SV^ = dt^^ + v'VJi^ ~ i'SI.V, (17) 

(quite obviously, since Vi — QijV-' and VkQij = 0). 



3 REVISITING AND REVISING THE SINGLE FLUID 
PROBLEM 

In order to lay the foundation for our analysis of the superfluid 
problem, it is useful to revisit the analysis of an ordinary perfect 
fluid star. In doing this we want to stay as close as possible to 
the equations used to describe the two-fluid problem. We know 
from previous work that the relevant model to compare to is a 
barotropic perfect fluid, eg. since the g-modes are absent from the 
pulsation spectrum of a non-rotating superfluid model (Lee 1995; 
lAndersson & ComejboOll : (Comeilf2002i) . Furthermore, we prefer 
to work with the number density n, the fluid velocity Vi and the 
chemical potential fi rather than the pressure P. 

3.1 The perturbation equations 

From thermodynamic principles we know that, for a barotropic or- 
dinary fluid we have E — E{n), and 



dfi = -dP . 

n 

This allows us to write the standard Euler equation as 

{dt + v^\/j)vi + + -l-) = 



(18) 



(19) 



where jl — /i/ma as before. In addition we have the continuity 
equation 

dtn + Vi{nv') = (20) 
and the Poisson equation for the gravitational potential 



V $ = 4nmBGn . 



(21) 



We now want to perturb these equations. First of all, conser- 
vation of mass for the perturbations is readily expressed as 

An = -nViC — >Sn^ -Vi(nf) • (22) 
Consequently, the perturbed gravitational potential follows from 



V^5$ = 47rGmB(5n = -47rGmBVi«') 



(23) 



In order to perturb the Euler equations we first rewrite Eq. <19> 



{dt + £v)v, + \/, (fi + ^ - -V ] = . 



(24) 



This form is particularly useful since the Lagrangian variation com- 
mutes with the operator dt + £v, i.e. 

A{dt + £v)v, = {dt + £v)Avi . (25) 
Perturbing J24t we thus have 

{dt + £v)Av, + V^ (A/i + A* - i A(t;^)) = . (26) 

We want to rewrite this equation in terms of the displacement 
vector ^. After some algebra one finds that 

dU^ + 2v'V,dt^^ + {v'Vjf^^ + V,;5$ + V,Vj$ 

-(V>f )Vj/i + V,A/i = 0. (27) 

Finally, we need 

' dp,'^ 



A/i = + C V»A = (1^) Sn + C V, 



(2t 

\dn 



(28) 



Given this, we have arrived at the following form for the perturbed 
Euler equation 

dU^ + 2v'Vjdt^, + {v'VJf^^ + V,(5$ + C'V,Vj($ + fi) 



\ dn 



= . (29) 



This equation should be compared to Eq. (15) of 
iFriedman & SchutJ Il9783l . 



3.2 Conserved quantities: Tlie canonical energy/angular 
momentum 

Having derived the perturbed Euler equations, we want to construct 
conserved quantities that can be used to assess the stability of the 
system. To do this, we first multiply Eq. <29> by the number density 
n, and then write the result (schematically) as 



AdU + BdtS, + = . 



(30) 



We are omitting the indices, since there should be little risk of con- 
fusion. 

Defining the standard inner product 



(31) 



where the asterisk denotes complex conjugation, one can readily 
show that 



and 



{n,Bi) = ~{i,Br,T . (32) 



The latter requires the background relation Vi {nv^ ) = 0, and holds 
provided that n — > at the surface of the star. A slightly more 
involved calculation leads to 



In particular, we use 

' djl 



(33) 



dV = / V 



dy(34) 



which (again) holds as long as n — > at the surface, and 
1 



nrf'Vi5i,^dV 



imUBG 
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1 



'i-Km-g.G 



= / Vj ((5^$)*dl/ . (35) 



Inspired by the fact that the momentum conjugate to 5' is 
p{dt + f"* Vj)^', we now consider the symplectic structure 



r,,Adt^+-B^ 



(36) 



where rj and ^ both solve the perturbed Euler equation. Given this, it 
is straightforward to show that W{ri, ^) is conserved, i.e. dtW = 0. 
This leads us to define the canonical energy of the system as 

Ec = ^Widt^, = ^ m^, AdtO + (C, CO} ■ (37) 

After some manipulations, we arrive at the following explicit ex- 
pression 

+ (1^) M'^-^a^^f^dv (38) 

which can be compared to Eq. (45) o f'Friedman & Schutj 119783) . 

In the case of an axisymmetric system, eg. a rotating star, we 
can also define a canonical angular momentum as 



Jc 



Re {d^OAdt^+^B^ 



(39) 



The proof that this quantity is conserved relies on the fact that (i) 
VF(»7,^) is conserved for any two solutions to the perturbed Eu- 
ler equations, and (ii) dip commutes with pv^ \7j in axisymmetry, 
which means that if f solves the Eule r equations then so does d^O 
As elucidated bv lFriedman & Schutz U^78a b), the stability 
analysis is complicated by the presence of so-called "trivial" dis- 
placements. These trivials can be thought of as "integration con- 
stants" representing a relabeling of the physical fluid elements. A 
trivial displacement ^' leaves the physical quantities unchanged, 
i.e. is such that Sn — Sv^ = 0. This means that we must have 

Vi(pC) = 0, (40) 

{dt + £v)C = 0. (41) 

The solution to the first of these equations can be written 

pC = e'^"V,-Xfc (42) 

where, in order to satisfy the second equations, the vector Xk must 
have time-dependence such that 



{dt + £v)xk = , 



(43) 



This means that the trivial displacement w ill remain constant 
along th e background fluid trajectories. Or, as lFriedman & Schutd 
Jl978iJ) put it, the "initial relabeling is carried along with the un- 
perturbed motion." 

The trivials may cause trouble because they affect the canon- 
ical energy. Before one can use the canonical energy to assess the 
stability of a rotating configuration one must deal with this "gauge 
problem." The way to do this is to ensure that the displacement 
vector ^ is o rthogonal to all trivials. A pr escription for doing this is 
provided bv 'F riedman & Schutd ^978?). In particular, they show 
that the required canonical perturbations preserve the vorticity of 
the individual fluid elements. Most importantly, one can also prove 
that a normal mode solution is orthogonal to the trivials. Thus, nor- 
mal mode solutions can serve as canonical initial data, and be used 
to assess stability. 



3.3 Example: Instabilities of rotating perfect fluid stars 

The importance of the canonical energy stems from the fact that it 
can be used to test the stability of the system. In particular, we note 
that: 

• Dynamical instabilities are only possible for motions such that 
Ec = 0. This makes intuitive sense since the amplitude of a mode 
for which Ec vanishes can grow without bounds and still obey the 
conservation laws. 

• If the system is coupled to radiation (eg. gravitational waves) 
which carries positive energy away from the system (which should 
be taken to mean that dtEc < 0) then any initial data for which 
Ec < will lead to an unstable evolution. 

Consider a real frequency normal-mode solution to the per- 
turbation equations, a solution of form ^ — ^g»("'+'"^) Q^e can 
readily show that the associated canonical energy becomes 



Ec 



io{OAO~,:{OBO 



(44) 



where the expression in the bracket is real valued. For the canonical 
angular momentum we get 



Jc 



(45) 



Combining Eq. j44t and Eq. j45t we see that, for real fre- 
quency modes we will have 



Ec — Jc — OpJc 

m 

where Op is the pattern speed of the mode. 

Now notice that Eq. <45t can be rewritten as 

Jc ^ {^ipv-VO 

— —mu) + m- 



Using cylindrical coordinates, and 



(46) 



(47) 



Q.ip^ , one can show that 



But 

lit xi).\ ^ icT 

and we must have (for uniform rotation) 
1 



1 + 



Jc/m^ 

^ — — ^ (Tb 



n 1 



1 

m 



(48) 



(49) 



(50) 



Eq. <50> forms an integral part of the 'Friedman & Schut3 
( 1978b) proof that rotating perfect fluid stars are generically unsta- 
ble in the presence of radiation. The argument is as follows: Con- 
sider modes with finite frequency in the limit. Then Eq. <50> 
implies that co-rotating modes (with Op > 0) must have Jc > 0, 
while counter-rotating modes (for which ap < 0) will have Jc < 0. 
In both cases Ec > 0, which means that both classes of modes are 
stable. Now consider a small region near a point where ap = (at a 
finite rotation rate). Typically, this corresponds to a point where the 
initially counter-rotating mode becomes co-rotating. In this region 
Jc < 0. However, Ec will change sign at the point where ap (or, 
equivalently, the frequency uj) vanishes. Since the mode was stable 
in the non-rotating limit this change of sign indicates the onset of 
instability at a critical rate of rotation. 

3.4 Example: Tlie r-mode instability 

In order to further demonstrate the usefulness of the canonical en- 
ergy, let us prove the instability of the single-fluid r-modes. 
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F or a general inertial mode we have rcf. lLockitch & FriedmanI 
il999l) who provide a discussion of the single fluid problem using 
notation which closely resembles the one we adopt here) 



and 



5^ ^5n^ Q,^ 



(51) 



If we also assume axial-led modes, like the r-modes, then we have 
5vr ~ and the continuity equation leads to 

V ■ 5i7 ~ ^ V ■ C ~ • (52) 
Under these assumptions we find that Ec becomes (to order 



p [\dte - \v ■ VCI' + r*^' V,V,(<E. + p)] dV .(53) 



We can rewrite the last term using the equation governing the ax- 
isymmetric equilibrium. Keeping only terms of order we have 



4 THE SUPERFLUID PROBLEM 

In this Section we generalise the Lagrangian perturbation analysis 
to the two-fluid model for superfluid neutron stars. In order to sim- 
plify matters we only consider the case of vanishing entrainment. 
The, significantly more complicated, general case will be discussed 
in a subsequent paper. 



4.1 The perturbation equations 

As in the single-fluid problem, we begin by deriving the equations 
governing Lagrangian perturbations of the system. Assuming van- 
ishing entrainment, i.e. letting a = 0, we have the Euler equations 

1 2' 

where we recall that x = n or p. Clearly, we must introduce 



(64) 



(54) 



two distinct Lagrangian displacement vectors ^x- To distinguish be- 
tween the two possibilities we use variations Ax such that 



A bit more work then leads to 



AxQ = SQ+£i,^Q . 



(65) 



-!:!"f*C'V>V,(r sin^ 
and 



in'] (55) 



The analogous generalisation of the single-fluid formulae for co- 
and contravariant vectors is straightforward. 

The two Lagrangian variations are naturally introduced in 
such a way that 



Iw'V^Cjf =f^' {m^l^l' -2imr' sin0cos6l [C^^* - C^'*] 



(66) 



+ [cos^ + sin^ l'] } (56) which leads to the perturbed continuity equations taking the form 



which means that the canonical energy can be written in the form 



AxTix = -nxVifi — > 5nx = — Vi(nxCx) 



(67) 



Ec~--J p{{mn~uj){mn + uj)\^f 

-2imQ.'^r^ sin 6 cos 6 [C^^"^* - ^"^5**] } dV 

for an axial-led mode. 

Introducing the axial stream function, 

iU 



where Yr = 
have 



iU 



d\rm iujt 



\OeYi e 



(57) 

(58) 
(59) 



With these definitions, it is very easy to derive the perturbed 
Euler equations. Simply comparing Eq. <24t to Eq. {(A\ we see that 
we must have 

-(V«Cx)Vi/Ix + V,Ax/ix = 0. (68) 
To express this in terms of the displacement vectors we need 

Ax/ix = <5/ix + CiVi/ix 



( 9wx ) „^ 



Yl^{9, ip) are the standard spherical harmonics, we 
1 



+ 



sm 



SUy + ^iVi/ix 



(60) 



and 

ir' sin^cose* - C'^C"*] 

1 cos 9 



and we arrive at the following form for the perturbed Euler equa- 
tions 



+ 2viV,dtC + {viVjfQ + Vi5$ + C^V,V,($ + Ax) 



smf 



m\U\'[YrdeYr* + YrdeYr 



(61) 



(I:) 



0. (7( 



After performing the angular integrals, we find that 



1(1 + 1) 



1(1 + 1) 

Comb ining this with the r-mode frequency iLockitch & FriedmanI 
Il999h 



From this equation it is clear that the two fluids are coupled. 
2 In order to proceed we need to understand the nature of this cou- 

p\U\ dr (62)pjjj^g better. In particular, we note that the perturbed gravitational 
potential depends on both displacement vectors. We have 



Lu — mO, 



1 



1(1 + 1) 



V^<5-I' = 47rmBG(<5n„+5rip) = -47rmBGV.(nnCn+npC) -(71) 
Since this is a linear equation we can write the solution as 



(63) 



we see that E^ < G for alU > 1 r-modes, i.e. they are all unstable. 
The I = m = 1 r-mode is a special case, leading to Ec — 0. 



(5$ = + S^p 

where we define 



x=n,p 



(72) 



6 A^. Andersson, G.L. Comer and K. Grosart 



V (5<l>x = 47rmBG(5rix 



-47rmBGV.(nxCi) 



(73) 



In analogy with the single fluid case, we can write the per- 
turbed Euler equations in the schematic form (after multiplying 
Eq. Go} by n^) 



Ax9t + BxatCx + CxCx + DxCy = 



(74) 



It should be noted that the first three terms are obvious generalisa- 
tions of the single fluid case. Now they pertain to each of the two 
fluids. The last term is new, and describes the coupling between the 
fluids. Explicitly, it takes the form 



-nxVi 



9Mx 



H-nxVi^^y 



(75) 



From this we see that the fluids are coupled (i) "chemically" 
through the equation of state, and (ii) "gravitationally" because of 
the fact that variations in one of the number densities affects the 
gravitational potential, which then influences the other fluid. 

4.2 Conserved quantities in the superfluid case 

We want to derive conserved quantities similar to those in the 
single-fluid case. To do this we again use the inner product. Given 
the results from Sec.|3| it is easy to show that we have the following 
symmetries 



(77x,^xi^x) = {Cx,Axr?x)* , 
{»?x,-BxCx) = -(Cx,-Bx??x>' , 
(??x,CxCx) = (Cx,Cx??x>* , 



(76) 
(77) 
(78) 



where rj^ can (at this point) be any vector field. 

Next we want to introduce symplectic structures that would be 
natural generalisations of the one we used to construct the canoni- 
cal energy and angular momentum in the single fluid problem. To 
do this we consider two sets of solutions [Cm Cp] ™d ['yn, ?7p] to our 
perturbation equations. Then we define 



(79) 



Given this definition and the above symmetry relations, it is 
straightforward to show that 



Analogously we introduce 



(80) 



W^p('7p,Cp) 



ApdtVp + i^Bprip,ip 



which leads to 

dtWp = -{71p, Dp^n) + (L>p7,„, 7^ 



(81) 



(82) 



Intuitively one would expect the sum Wn + Wp to be con- 
served. That is, the coupling terms in Eq. <70> should facilitate 
non-dissipative energy transfer between the two fluids. We will now 
prove that this is, indeed, the case. 

Explicitly we have 



dtWn 



nn77n* Vi 



dV 



dfln 

dup 



V]{npr)l* 



dV 



nnr]"V^5i_^^pdV + / nnCnV,(5^p'l>p)'dl/ (83) 



and 

dtWp 



UpTlpVt 

jipCpV 



Up 

dfLp 



dV 



dV 



- J npTjl*V,5^^^ndV + J npC;V,(5^„$n)*dl^ . 
Consider the first two terms of <84> . Using the fact that 

/9Ap\ 
V dun J „„ 



dup 



and assuming that Un and Up both vanish at the surface^ 



(84) 



(85) 



i* 5-7 

npTip V, 



dV 



and 



'^pCpV^ 



djli, 
drir: 



Vj {ripril* 



dV 



V,(npC^) 



dV 



(86) 



dV 



(87) 



(88) 



This means that, when the two expressions Eq. <83> and Eq. <84t 

are added, the first two terms of each expression will cancel each 
other 

To rewrite the terms involving the gravitational potentials in 
Eq. J83> and Eq. J84> . we need to use 

J nnVpV^{5^p^p)dV = J npC;V.(5^p$n)*dV 

Then it is easy to show that the last two terms in Eq. <84t will cancel 
the last two terms in Eq. <83> once we add the two expressions. 

In other words, we have proved that, in addition to <76t - <78> . 
the superfluid equations have the following symmetry 

{Vn, Dn^p) + (77p, DpCn> = (Cn, D^Vp)* + (Cp , DpT^r.)* (89) 

and it follows that 

W = Wn{Vn,^n)+Wp{r^p,^p) (90) 

is a conserved quantity. 

By analogy with the single fluid case, it now makes sense to 
define the canonical energy of the system as 

me 



[Wn{dt^n,^n) + Wp{dt^p,^p)] 



(91) 



^ Somewhat artificially, we assume that the rotating background is such 
that the two fluids have a common surface. In reality, the outer layers of a 
neuti'on stai' will not be superfluid and one would have to add a single fluid 
envelope to our model. The analysis of such composite models is beyond 
the scope of the present analysis. 
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This can be written 



+ (Cn,C„?„>* + (Cp,CpCp)* + (Cn.^nCp)* + (Cp.I^pCn)*} ■ (92) 

After some manipulations, we arrive at the following final explicit 
form 

[P^dtinf + PpI^CpI' - PuKV,Q\^ 



KdUn 



ISnJ 



dUr, 



\Snp\ 



9nc 



It is worth noting that the symmetries of our system of equa- 
tions imply the existence of a quadratic Lagrangian for the pertur- 
bations, which in turn implies the conservation of 14^ = Wn + Wp. 
We can derive the superfluid equations of motion i74i from a vari- 
ational principle whose action is (for clarity assuming that the dis- 
placement vectors are real, the generalisation to the complex case 
is straightforward) 

/ = J £dV^^ [{in,Anin) + (Cp, Ap?p) + (L BnCn) 

+ (Cp)-BpCp) - (CniCnCn) ^ (Cp,CpCp) 

-(ep,A,Cn)-(?n,DpU] • (94) 

Then the momentum conjugate to follows from dC/d^x, while 
the equations of motion <74t can be derived from the standard 
Euler-Lagrange equations 



dt 



dC 



dC 



. 



Finally, the canonical energy for the system is 



■i dC -i dC 
f ^ V 



c 



dV 



(95) 



(96) 



One can readily verify that these formulas lead to the given results. 

In an axisymmetric system we can also define a conserved 
angular momentum; 

Jc = -^w.(a^Cn,Cn) - ^Wp(a^Cp,Cp) 



-niB Re ^(^d^^n, Andt^n + ^-BnCn^ 

pCp,^p9tCp + iBpCp)} • (97) 



Now one can readily use Eq. <92> and Eq. <97> to prove that, 
for a normal mode solution to the problem, = ^,^e*''"'^"''"^*\ the 
canonical energy and angular momentum will still be related by 
Eq. <46t . Furthermore, it follows that we must have Ec — Jc = 
for dynamically unstable (complex frequency) modes, just like in 
the single fluid case. This is hardly surprising, but it could turn out 
to be a very useful result. One could, for example, hope to be able 
to use our expressions for the canonical energy and angular mo- 
mentum t o derive necessa ry criteria for the superfluid two-stream 
instab ility iAndersson et al. 2003 ; Prix et al. 2004 ; Andersson et al. 
l2002h. 



4.3 Trivial displacements 

Before we close this section, let us address the issue of trivial dis- 
placements in the two-fluid problem. We clearly need to deal with 
two sets of trivial displacements, one for each fluid degree of free- 
dom. Fortunately, the analysis of these trivials is essentially iden- 
tic al to that of the single barotr opic fluid case discussed in Sec.|3| 
see iFriedman & SchutzHl9783) for further details. 

In the superfluid problem, the trivial displacements are such 

that 



5nx ~ 5vf — . 



(98) 



It is easy to see from the superfluid perturbation equations that the 
two sets of equations that determine the functional form of the triv- 
ials are identical to the corresponding single fluid equations. This 
mean^StBi^t the single fluid result can be adapted to the superfluid 
problem: we simply have one trivial displacement ((^x) per fluid. 

In our case, we can take the requirement that the canonical 
displacements ^x must be "orthogonal" to the trivials to mean that 
we should have 



H/„(Cn,Cn) = vyp(Cp,ep) = o 



(99) 



This condition leads to the trivial displacements having identical 
form to those of the single-fluid problem. However, it is worth 
noticing that we could in principle permit the somewhat less re- 
strictive condition 



W„(Cn,Cn)+Wp(Cp,ep)=0 



(100) 



and still ensure that the trivial displacements do not affect our con- 
served quantities. We have not yet investigated the implications of 
this possibility. 

Final ly, and most importantly, one can readily extend the cal- 
culation of lFriedman & SchutJ il978ah to prove that normal modes 
are (usually) orthogonal to the trivials also in the present case. In 
other words, normal modes may serve as canonical initial data. This 
is extremely useful as we hope to use the canonical energy and an- 
gular momentum to assess the stability of various superfluid normal 
mode solutions. 



5 INSTABILITIES OF ROTATING SUPERFLUID STARS 

The main motivation for the present investigation was the lack of 
proper instability criteria for rotating superfluid stars. The need 
for such criteria is clear given that it has long been acknowl- 
edged that the astrophysical relevance of the gravitational-wave 
driven CFS instability (of both f- and r-modes) may depend on 
the extent to which superfl uid dissipation [like mutual friction 
iLindblom & Mendelf 1995, 2000)] counteracts the growth of the 
unstable mode. It would seem obvious that, before worrying about 
such issues, one ought to establish that the instability is actually 
present once the star becomes superfluid. We are not aware of any 
such proof, despite the number of investigations of unstable oscilla- 
tions of superfluid stars that exist in the literature. It should be clear 
that this is a non-trivial issue given the simple fact that the two 
fluids are only weakly coupled, and may in fact rotate at different 
rates. 

In this section we provide the first applications of our con- 
served canonical energy and angular momentum for the two fluid 
problem. We first present an argument in favour of the simple single 
fluid criterion for CFS instability — that the instability sets in when 
the pattern speed of an originally backwards moving mode passes 
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through zero in the inertial frame — holding also in the superfluid 
problem. Having done this, we discuss the superfluid r-mode insta- 
bility. 



This concludes the proof that the criterion for the onset of radiation 
driven instabilities of modes that have a finite frequency limit as 
f2 — > remains as in the single fluid case. Modes become unstable 
when the (inertial frame) pattern speed changes sign. 



5.1 The superfluid CFS instability 

The main question here is: Does the simple criterion that a counter- 
rotating mode becomes unstable when the pattern speed changes 
sign remain valid also in the two-fluid problem? Intuitively, one 
might expect this to be the case, but it nevertheless warrants a proof. 

We approach the problem as in the single fluid case. Assum- 
ing a real frequency mode solution to the perturbation equations, 
Eq. 197} leads to 



Jc = -m {i^[(Cn, PnCn) + (Cp, PpCp)] 

- [(Cn, ipuVu ■ V^n) + (^p, ippVp ■ V^p)]} 

Divide through to get 

Jc 



(101) 



{^n,PnCn> + (?P,PpCp> 



-mto + m 



{^n,ipnVn ' Vgn) + {^p,iPpVp ' V^p) 



(102) 



Using cylindrical coordinates, we see from the results in 
Sec.l33lthat 



■ iPnCiviVjCn = Pn^nlml^nf + i{^n X (,n) z] 



(103) 



6 THE SUPERFLUID R-MODE INSTABILITY 

To conclude this paper we will discuss some aspects of the r-mode 
instability for rotating superfluid stars. Although we are not yet at 
a point where we can discuss the general problem (since we did 
not include entrainment in our derivation), we can still learn quite 
a lot about the issues that arise when we consider two coupled flu- 
ids. Furthermore, we are not aware of any previous proof of the 
presence of an instability in the case when the two fluids rotate at 
different rates. The nature of the various inertial modes of oscilla- 
tion (of which the r-modes form a sub-class) of a superfluid star 
has, however, been discussed by several authors. We will draw on 
these investigations for information concerning the nature of the r- 
modes in both (i) the general case of a background star such that 
the two fluids rotate at different rates (although with respect to the 
same axis), and (ii) the special case of co-rotating fluids. 



We then know that 

Pn!^n(m-l)|Cn|' < -tpnCm^'^^jC < /On^n (m + 1) j^n | ^ ( 104) 

and it is easy to see that we will have^ 

(Cx, ipxWx • VCx) > -(m + l)t^xax (105) 
where we have defined — (^x, Px^x) > 0. We also get 
(Cx, ipxWx • V^x) ^ -(m - l)f2xax ■ (106) 
These results can be summarised as 

crp{an + ap) — ^1 + — ^ [Slnan + Slpfflp] ^ Jc/m^ 

^ (Tp{an + ap) - f 1 - —1 [finfln + ^^pOp] .(107) 
Z \ 771/ 

This is the relation we need. Provided that Slnin + SlptSp > 
(which is obviously true if the fluids both rotate in the positive di- 
rection) we easily show that 

• if we let [fin, ^p] — > while lo is finite, then for co-rotating 
modes we have ap > 0, which means that Jc > 0. In contrast, for 
counter-rotating modes ap < and Jc < 0. In both cases, we will 
get Ec > which indicates that all (finite frequency) modes are 
stable. 

• if we consider a region near = for finite rotation rates, 
then ap = implies that Jc < 0. This means that the mode is 
stable as long as ap < 0, but when ap changes sign (and the mode 
becomes co-rotating) we will have iSc < and an instability. 



We are assuming that the two fluids rotate in the same direction here. Al- 
though it is possible to construct models for counter-rotating backgrounds 
these cases are somewhat pathological, and we do not expect them to have 
any astrophysical relevance. 



6.1 The case with relative rotation 

The g eneral case, in which Qn ^p, was recently discussed by 
IPrix et al. 1 2004) . From that study we learn that, in absence of en- 
trainment, the r-mode fluid motion must be such that only one of 
the fluids oscillates. This means that we will have two classes of 
modes, corresponding to 



(5«x / , SVy 



. 



mf^x 



1(1 + 1) 



(108) 



Since this implies that only one of the two displacements vec- 
tors is non- vanishing, it is straightforward to show that the canoni- 
cal energy Eq. <93t reduces to 

= i y {pxi^tCxi^ - Pxi«^ Vj,tf 1^ 

+pxa?x*V,V, ($-f /ix)}dV (109) 

where we have assumed that Sn^ and 5$ are of higher order in 
the slow-rotation scheme (as in the single fluid case), and also used 
= w*. Noticing the close resemblance of this result to the expres- 
sion for the case of a barotropic single fluid, Eq. <53L we readily 
infer that both these classes of modes will be unstable due to the 
emission of gravitational radiation. 

However, this result is likely of very limited relevance. Even 
though it is expected that the two rotation rates will be slightly 
different in astrophysical neutron stars, the entrainment coupling 
will affect the r-modes significantly [as elucidated bv IPrix et alJ 
(2004)], naturally leading to both displacements being non-zero. 
This means that, despite being of conceptual interest, the case we 
have discussed here is pathological. 

6.2 The co-rotating case 

In the special case of = Op = Q, we know 

( Andersson & Comet 20o H: IComed l200l Lee & Yo shida 20^ 
lYoshida & Lea .2003b: .Prix et alj|2004 l) that there will exist two 
classes of r-modes. One is such that the two fluids move in phase, 
while the other has the two fluids counter-moving. Furthermore, for 
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non-stratified stars one can show that these two degrees of freedom 
decouple. Given the relative simplicity of this case, we will focus 
our attention on it. The general case essentially follows as a linear 
combination of the two results we present. 

We begin by introducing the two classes of displacements 

= !^^^l±Ih^ (110) 

Tin + Up 

and 

er = e"-er- (hd 

These are clearly such that, when ^" = and the two fluids move 
together only is present, while the total momentum flux van- 
ishes when — 0. 

Now consider purely axial r-mode solutions such that ^~ — 0. 
In this situation the canonical energy Eq. |93) can be written as 
(since vi = = v-') 

Ec^\j {p\dte? - pWV jit? + Pi\ii*y^y 

K V> V.Mn + np V> Vj^ip] ]dV (112) 

where we (again) neglect the higher order contributions from 5n^ 
and 

In this expression, the last term can be rewritten using the fact 
that we must have 

Vi/in = Vi/ip = Vi^i (113) 

if the two fluids rotate at the same rate, cf Eq. <64t . This immedi- 
ately leads to 



where Xp = ^ is the proton fraction. Since the expression in the 
bracket has the same form as in the single fluid case, and the pref- 
actor pnXp is positive definite, it is easy to prove that Ec < Q also 
for these counter-moving modes. 

6.3 Gravitational-wave emission 

At this point it is appropriate to discuss how efficient the oscilla- 
tions of a superfluid star are as a source of gravitational waves. Af- 
ter all, we have shown that the superfluid r-modes generally lead to 
Ec < 0, and should therefore be driven unstable by radiation. The 
relevance of the instability is then largely dependent on the rate at 
which the motion generates gravitational radiation. 

To address this problem, we consider a source with weak inter- 
nal gravity, and focus our attention on a single pulsation mode with 
time-dependence exp(iLjt). We also assume that the background is 
such that the two fluids rotate at different rates: and respec- 
tively. 

T he gravitational-wave luminosity follows from e.g. lThorn3 
h980t). 



(120) 



dt 
where 

Ni = 



1=2 



4nG {l + l)(l + 2) 



(121) 



c2'+i/(/-l)[(2/ + l)!!]2 ■ 

The first term in the bracket of Eq. <I20t represents radiation due 
to the mass multipoles. These are, quite generally, determined from 



(where n — Un + rip) and 



(114) 



SDi^ = / TooYil.r'dV 



Ec 



P{\dti 



+ |2 



'v,c+l'+CUl*v,v,($ + A)}dy 



Clearly, this result is (provided that we identify f with identical 
to that of the single fluid problem, Eq. <53> . Hence, the instability 
of pure modes follows from the calculation in Section 3.4. This 
is not very surprising given that the degree of freedom we are con- 
sidering is such that the two fluids move together. 

Next we consider the canonical energy for counter-rotating 
modes, which are such that 



(122) 



where T^j^ is the contribution to stress-energy tensor associated 
.(115)ivith non-axisymmetric motion in the source. The second term in 
the bracket of Eq. <I20t corresponds to the current multipoles, 
which follow from 



rinC + ripCp = . 

In this situation only ^~ is present, and we have 

,11 _ Tip ^_ 
Si Si 5 



Si 



' _ Si _ 



(116) 

(117) 
(118) 



(123) 



where oc f x VY/" are the magnetic multipoles iThomd 
.1980 ). 

By taking the Newt onian limit of the relativistic stress-energy 
tensor (see, for example. IComed i2002h ) 



Til' ,T, rJ^ I I I' 

M = *<5m + P Xi^+n Pu 



(124) 



Incidentally, the latter of these relations emphasizes the fact that 
pure ^~ modes can only exist for non-stratified stars. We know 
from the results of Prix et al. ( 2004) that ^" and will have the 
same functional dependence on the radial coordinate r. This means 
that, in order for both displacements to be non-zero, they must be 
proportional. This is clearly only possible if n-a/up is constant. 

Using the above relations in the expression for the canonical 
energy Eq. <93t . one can show that 



[which is easily done usin g formulas given in Appendix A of 
lAndersson & Comej|200ll) 1. one can show that 

Too~{p^ + Pp) (125) 
and 

Toj « p„Vn + /9p«p (126) 

for a Newtonian source. Perturbing these expressions we find that 



5Dlm = 

and 

5 Jim ~ 



{Spn + Sppy'Vii^dv 



(127) 



l + l 



r {pnSvn + ppSvp 



(119) 



+ 5pn^n + SppVlp) ■ YiJ'dV 



(128) 
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As discussed in the previous section, it is sometimes instruc- 
tive to express the superfluid formulas in terms of the variables 



5^ = P^5v^ + P^s^ 
P P 

and 

Using these, together with 

5p = 5pp + (5pn 

we get 

SDi^ = j Spr'Yi*^dV . 
and 



(129) 



(130) 



(131) 



(132) 



SJlm — - 



l + l 



+5p„{Qn-np)]-YiZ*dV 



(133) 



When written in this form, the formulas closely resemble the stan- 
dard single-fluid results. The only new feature is the last term in the 
bracket of Eq. <133t . 

An interesting question concerns whether it is possible to have 
oscillations in a superfluid star that (at least at this post-Newtonian 
level) do not radiate gravitationally. From the above equations we 
immediately deduce that in the case of a non-rotating star (with 
f2n = f^p = 0) or a co-rotating star (when — f2p) we must 
have 5p — 5v^ = in order not to have any gravitational radiation 
emission. This result is quite intuitive since, as discussed in the pre- 
vious section, the co-moving degree of freedom represents the total 
momentum flux. It is not surprising to find that motion which cor- 
responds to zero momentum flux does not radiate gravitationally. 

However, when combined with the results obtained from our 
superfluid canonical energy, these results illustrate some of the 
complexities associated with a discussion of gravitational-wave 
driven instabilities in a superfluid star. In particular, we con- 
clude that even though they are formally unstable (Ec < 0) any 
purely counter-moving modes (for which — 0) will not grow 
(dE/dt = 0). 



7 CONCLUDING REMARKS 

With this paper we have taken the first steps towards a Lagrangian 
perturbation framework for rotating non-relativistic superfluids. 
The primary motivation for this work, and the key application con- 
cerns the stability properties of rotating supe rfluid neutron stars. 
Our analysis generalises the classic work of iFriedman & Schutd 
iT978a.b) to the case of stars which require a multi-fluid descrip- 
tion. We have applied our framework to the problem of dynami- 
cal and secular instabilities of a simplified superfluid neutron star 
model in which the two fluids are allowed to have different uniform 
rotation rates around the same axis, and where the entrainment ef- 
fect is neglected. We have demonstrated that the criterion for the 
onset of radiation driven instabilities for modes that have a finite 
frequency as the background rotation vanishes remains unchanged 
from the ordinary fluid case. We have also considered the super- 
fluid analogue of the r-mode instability, and found that both co- and 
counter-moving superfluid r-modes have negative canonical energy 
and will therefore be driven unstable by gravitational-wave emis- 
sion. 



Of course, our neglect of the entrainment is a serious limita- 
tion of the formalism. However, as we stated earlier, if we could 
not make progress in the case of vanishing entrainment it would 
be pointless to proceed to the general case. Naturally, we are en- 
couraged by the progress we have made, which provides crucial 
benchmarks for the onset of instability in multi-fluid systems. In 
the near future, we will aim to generalize our results to include en- 
trainment. This is absolutely essential for the study of gravitational- 
wave driven instabilities, since the primary superfluid damping 
mechanism is the mutual friction which depends crucially on the 
entrainment. The inclusion of entrainment will also be important 
for future studies of the superfluid two-stream instability, since en- 
trainment may be the dominant coupling to drive the instability in 
the interiors of neutron stars. 
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